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Abstract -A modified spin-fermion model is proposed to describe the physics of the underdoped 
phase of the t — J model of strongly correlated electrons. 


More than a decade ago, a spin-fermion (SF) model was 
put forward in an attempt to account for the observed 
non-Fermi liquid anomalies in the doped cuprates near 
optimal doping. [1] It describes low-energy fermions with 
a large Fermi surface (FS) interacting with each other via 
soft collective spin excitations. The bare spin propagator 
emerges from the integration of the high-energy fermionic 
modes, 

xo{q,uj,,)^{‘4 + {q-Qf + GF)-\ ( 1 ) 

where a;„ = ^'nnjfi are Matsubara bosonic frequencies. 
The antiferromagnetic (AF) spin correlation length, 
and the spin velocity, Vg, are the high-energy phenomeno¬ 
logical parameters. The spin velocity is of order of the 
Fermi velocity, vp. 

At criticality (^af = oo), the singular nature of the low- 
energy AF paramagnons singles out the so-called hotspots 
on the FS. These are the points on the FS for which the AF 
momentum Q = Qaf = (F, 7 r) coincides with 2kF- Due 
to the enormous enhancement of the AF fluctuations at 
or near the quantum critical point (QCP), the fermions in 
the vicinity of hotspots interact strongly with each other 
via singular AF paramagnons. This in turn results in a 
strong damping of the propagating spin modes in the low- 
energy limit. The spin modes become dissipative and the 
fermion quasiparticles acquire a short lifetime. As a result, 
the electron system displays a singular non-FL behavior 
at low energies. m 

There are, however, two necessary conditions for this 
approach to work. First, a SF model must be at weak 
coupling. As a result, a spin-fermion coupling g should 
be smaller than the fermionic bandwidth. Only in this 


case can one perturbatively integrate out high energies to 
obtain an effective low-energy action in which the spin 
degrees of freedom are represented by the bare spin prop¬ 
agator O- Second, the low-energy one-loop spin-fermion 
dynamics renormalizes the bare spin propagator produc¬ 
ing the effective susceptibility, 

Xeff{q,i^n)'-^{—^ + iq-Q)‘^+fAF+T\^n\) ( 2 ) 

where a dissipation coefficient j ^ g'^. The non-FL low- 

energy behaviour sets in provided 7 is large enough to 

2 

allow for the dropping of the ’’ballistic” piece ^ in the 
presence of a much more pronounced ’’dissipative” con¬ 
tribution. These two conditions formally contradict each 
other, however. 

A possible way out of that is to employ a 1/A^ expan¬ 
sion, where N stands for a number of the fermion flavours 
with A^ = 1 being the physical case. It can be shown that 
a random phase approximation (RPA) generated by a one- 
loop particle-hole expansion contributes to the order 0 ( 1 ) 
whereas all other higher loop diagrams vanish in the limit 
N -F 00 . At the same time, 7 remains large and finite. 
The 1/A is then the parameter that keeps the theory un¬ 
der control. The SF theory was developed in detail, along 
these lines, in a seemingly self-consistent way in Ref. [T]. 
A good agreement with experiment was declared as well. 

However, such a theory is incomplete. The 1/A ap¬ 
proach can only be applied reliably, provided the physics 
remains qualitatively the same at any A. This is not the 
case in the present situation. Moreover, the 1/A approxi¬ 
mation fails to take proper account of the non-FL physics 
displayed by the underdoped cuprates. [3] More specif¬ 
ically, it was recently shown that the 1/A approach to 
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a 2 + 1 SF model misses some important contribution in 
both perturbation and renormalization group theory 113, 
which invalidates its application in this context. The pro¬ 
posed 1/A^ SF model is in this way not well defined and 
this makes a direct comparison with experiment incom¬ 
plete and unjustified. 

A new interesting approach to treat a strong Landau 
damping within a framework of a weakly coupled SF 
model in a controlled way was recently proposed in [3. In¬ 
stead of the introduction of the N fermion flavours, the au¬ 
thors consider a reduced FS to account for the low-energy 
dynamics of the optimally doped cuprates. Such a FS 
comprises eight hotspots which are essentially the inter¬ 
sects of the non-interacting hole-like FS with the magnetic 
zone boundary. [3 This zone boundary is fully sensitive to 
umklapp scattering processes which gap out the FS com¬ 
pletely at half-filling. A small angle S between the Fermi 
velocities of two hot spots involved in scattering processes 
is introduced in this new approach. The Landau damp¬ 
ing modifies the form of the bare spin propagator as given 
by (l2|) with 7 being now given by 7 g^jsinS. In the 
limit (5 <C 1, the Landau damping is strong even within a 
’’weakly coupled” SF theory. This keeps the theory under 
control and it is, in principle, in agreement with experi¬ 
mental data. In particular, the proposed theory features 
an instability towards a formation of a quadropole charge 
order (CO). The magnetic-mediated interaction is shown 
to give rise to charge bond order with the momentum di¬ 
rected along the Brillouin zone diagonal. While the struc¬ 
ture factor is compatible with experiments, the magnitude 
and direction of the ordering momentum are not. In ex¬ 
periments, the charge ordering momentum appears to be 
always directed along either the horizontal or the vertical 
axis. 

By construction, the approach advocated in [3 relies on 
a specific shape of a reduced hotspot FS which should de¬ 
scribe the cuprates near optimal doping. Apart from the 
already mentioned stability problems of such a FS, it is 
not clear if this approach is applicable at lower dopings 
in which strong electron-electron interactions are known 
to be at work. The heavily underdoped region is char¬ 
acterized by strong electron correlations which keep the 
majority of the electrons well localized. Accordingly, the 
FS consists of small separated patches (electron pockets). 
It is still unclear how far a standard weak-coupling SF 
model can go to describe realistically the physics in this 
strongly correlated region as well. 

The goal of the present paper is to explore this issue 
further. We show that, at least close to half-filling, the 
standard SF model should be drastically modified to take 
into account strong electron correlations ~ a hallmark of 
the physics of the underdoped cuprates. 

The SF model can be derived formally [3 from the Hub¬ 
bard model which describes fermions hopping on a lat¬ 
tice with a tunneling amplitude tij and are subjected to a 


short-range Coulomb repulsive interaction U: 

Htu = (3) 

ka ^ 

Operator is a creation operator for electrons with spins 
a on site i and riia = cl^Cia is the on-site electron number 
operator with the spin projection cr. The next basic step is 
the Hubbard-Stratonovich (HS) decoupling of the quartic 
fermion interaction. The SF model is then described by 
the Lagrangian [3 

^tg — ^ ^ 

qa 

+ Y Xo\<l)$q$-q + 9 Y ^q-Q^-g- (4) 

<? Q 

where g ^ U is a, spin-fermion coupling constant. The 
electron spin operator 

^ ^ Cj^^Tg-g-' Cjg' , (5) 

(7,cr' 

couples to the auxiliary HS field ipi. In case the SF coupling 
5 ~ 17 is smaller than the fermionic bandwidth, LF ~ t, 
one can apply a standard perturbation theory controlled 
by a small parameter U/t. 

In the cuprates, however, the fermionic bandwidth is 
comparable to the Coulomb repulsion and this makes the 
weak-coupling theory problematic. Besides, close to half 
filling, the tunnelling amplitude gets effectively renormal¬ 
ized, i.e., t ^ teff xt, where a; -C 1 is the doping concen¬ 
tration. [3 This results in an extremely narrow bandwidth 
for the itinerant electrons which automatically become 
subject to the important (in this region) local no dou¬ 
ble electron occupancy (NDO) constraint. This makes the 
hopping of an electron on an already occupied lattice site 
extremely unfavorable since U 3> xt. This results in extra 
strong correlations between electrons in addition to those 
caused by the Pauli exclusion principle. The lightly doped 
cuprates are thus characterized by the strong electron cor¬ 
relations in the coupling regime U/teff ^ 1. Since the SF 
coupling g U, a straightforward extension of Eq.(|4]) in 
this regime automatically implies that g becomes large. 

Let us now see what happens in the limit of an infinitely 
strong Coulomb coupling. The Hubbard model goes over 
into 

H-U—ca = ^ (®) 

ija 

where the CiaS are the on-site Gutzwiller projected electorn 
operators: 

— TzGcrPz, Pi — 1 ('^) 

The model m exhibits the extreme case of strong corre¬ 
lations and is referred to as the U = oo Hubbard model. 
This model is certainly far from trivial and it admits an 
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exact solution only in Id. Under the assumption that, in 
the limit g oo, the SF model (|3|) is still meaningful, it 
must recover in this limit the physics displayed by Eq. (|6]). 
However, the limit g —>■ -l-oo pushes the energy levels of 
model to -|-oo and this destabilizes the system. We 
should therefore search for a modified SF model which is 
finite in the limit of strong electron correlations. 

To deal with the underdoped phase, it is more conve¬ 
nient to use as a reference framework the so-called t — J 
model of strongly correlated electrons: 

HtJ = ^ ^ ^ ^ ^ ^ SiSj, ( 8 ) 

ija ij 

where J = /U. The t — J model is essentially the 

strong coupling limit of the Hubbard model, i.e., Htjj = 
HtJ + o{J/t), t/U <C 1. It is believed to capture the basic 
physics of strongly correlated electrons in the underdoped 
cuprates. 

Within the recently proposed spin-dopon representation 
of the t — J model |1] , the local NDO constraint is enforced 
by the requirement of an infinitely large spin-dopon cou¬ 
pling. In this way, the t — J model turns out to be equiva¬ 
lent to a Kondo-Heisenberg lattice model of itinerant do- 
pons (doped holes) and localized lattice spins (localized 
electrons) at an infinite Kondo coupling: 

HtJ = Tijd\^dja + JSi ■ Sj + \ Si - Si, ( 9 ) 

ij(T ij i 

where T^- = 2tij + (3A/4 — g)Sij and the doer’s are the 
fermion operators to describe the dopons. The su{2) gen¬ 
erators Si describe the localized electron spins. The on¬ 
site itinerant dopon spin operator is given by Eq. where 
the electron operators being replaced by the dopon opera¬ 
tors. A global parameter A is to be sent to -|-cx) to ensure 
the selection of the appropriate physical subspace. The 
unphysical doubly occupied electron states are separated 
from the physical sector by an energy gap A. In the 
A —^ -too limit, i.e. in the limit that A is much larger than 
any other existing scale in the problem, those states are 
automatically excluded from the Hilbert space. 

In this limit, the high and low energy itinerant fermions 
cannot be separated out and this is another manifestation 
of the Mott physics. One cannot therefore integrate out 
high-energy fermions in that case. We assume however 
that this separation still holds in the spin-dopon represen¬ 
tation for well localized particles, i.e., the localized spin 
degrees of freedom. After all, this is believed to be the case 
for the undoped Mott insulator (the quantum Heisenberg 
model) - a 2d lattice of localized electrons. We believe that 
such a separation is valid at sufficiently low doping as well, 
since a large majority of the electrons remain localized in 
this regime. 

At half-filling, a: = 0, and model ([S]) reduces to the 2d 
AF Heisenberg model, 

j'^Si-Sj, J>0. 
ij 


This model exhibits long-range AF order at T = 0. How¬ 
ever small dopings {xc ~ 0.05) suffice to destroy that or¬ 
der leaving behind strong short-range AF fluctuations that 
persist up to optimal doping. This can be taken into ac¬ 
count as follows. We write 

Siir) = {-yn{xi,T)+rn{xi,T), (10) 

where fi(xi,T) is a smooth function of {x,t) with 
\n[xi,T)\ = 1. Here n describes the low energy AF fluc¬ 
tuations. The vector field m with \rn{xi,T)\ ^ 1 and 
n{xi,T) ■ m[xi,T) = 0 describes the ferromagnetic fluc¬ 
tuations which should in turn be integrated out. This 
approach is controlled by a l/ 2 s expansion, where s is the 
electron spin. Although the physical case corresponds to 
s = 1/2, the physics of the short-range 2d AF phase is 
correctly reproduced in this case. 

We then substitute the ansatz (ITO)) into ® and drop the 
small TO-dependent term in the spin-dopon interaction, 

^ ] Sj ■ Sj —j X 'y ^ • ■S-g- ( 11 ) 

i q 

The high-energy integration over m can then be readily 
done in a standard fashion m resulting in a SF La- 
grangian to describe strongly correlated electrons in the 
underdoped phase, 

Ltx = X! dqa{-iUJn + T;j{\))dq„ 

qa 

T ^ ^ Hq (9)''^5^-5 + a (12) 

9 9 

Here q = {q, iiOn) and 

3A 

r,-(A) = 2V+--^. (13) 

The spin propagator D[q) is given by Eq.®, where Vg ~ 
Ja is now the spin-wave velocity, and a is a lattice spacing. 

We thus arrive at the SF model at an infinitely strong 
coupling A. In contrast with the standard model ®, the 
renormalized hopping amplitude contains the spin-fermion 
coupling A as well, as given by Eq. dni). Because of this, 
in the physical limit, A -too, the model (jl2|) remains 
finite. In Id it reproduces the exact ground-state energy 
of Hamiltonian ®. m 

It is instructive to consider (HI at small values of 
X, X/t <C 1. This is clearly a unphysical assumption in 
the underdoped region and in the whole PG phase as well. 
In this hypothetical case, one can ignore the A-dependence 
of the hopping amplitude Ty. Equation (|12l) then reduces 
to the weak-coupling SF Lagrangian given by Eq.® with 
the identifications g ^ X and tj: 2t^. We thus see that 
ignoring strong electron correlations indeed results in the 
standard weak-coupling SF model. There is however a 
qualitative distinction between the standard representa¬ 
tion ® and that given by our equation (fT^ . The spin 
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velocity in Eq. m is of order of the Fermi velocity vp ^ ta, 
whereas the spin-wave velocity Vs in D{q) is determined 
by the dynamics of the localized spins, i.e., Vg ~ Ja. 

As a result, the correct energy scale in the problem is 
the spin fluctuation frequency, uig/ ^ Vs^~^/k, where k is 
an effective dimensionless spin-fermiom coupling which is 
essentially ^^/vg, and ^ := ^af- This frequency mea¬ 
sures the deviation from the QCP. At the QCP = oo), 
uigf = 0, and the system displays non-Fermi-liquid be¬ 
havior down to the quasiparticle fermionic energy w = 0. 
Away from the QCP, the correlation length ^ is finite. For 
small energies, w << ujgf, the system displays FL behav¬ 
ior, although with a damping term scaling inversely with 
uigf , rather than with the Fermi energy as in conventional 
metals. At frequencies above ojgf, the system crosses over 
into the non-Fermi liquid quantum-critical regime. In the 
conventional SF approach, ujgf ~ whereas in our 

case ujgf ~ Since Wgf{vg)/wgf{vF) = {J/tY « 1) 

this implies that the energy window within which the es¬ 
sential non-Fermi-liquid behavior sets in is in fact much 
larger than that predicted by the conventional SF ap¬ 
proach, solely based on the itinerant fermion description. 

To conclude, we propose an effective SF model that 
takes into account strong electron correlations. These cor¬ 
relations bring us naturally to the SF model at infinitely 
strong coupling with a simultaneously renormalized hop¬ 
ping amplitude. This theory is difficult to implement an¬ 
alytically, since there is still no machinery available with 
the proper technical tools. However, it is just this effective 
model that is appropriate to address the cuprates at very 
low doping, rather than its weak-coupling SF version. The 
latter leaves out the effects produced by the NDO - the 
essence of the physics of strongly correlated electrons. 
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